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Beam stability and halo formation in high-intensity axi- 
symmetric 2D beams in a uniform focusing channel are ana- 
lyzed using particle-in-cell simulations. The tune depression 
- mismatch space is explored for the uniform (KV) distribu- 
tion of the particle transverse-phase-space density, as well as 
for more realistic ones (in particular, the water-bag distribu- 
tion), to determine the stability limits and halo parameters. 
The numerical results show an agreement with predictions of 
the analytical model for halo formation 



I. INTRODUCTION 

There is an increasing interest in high-current applica- 
tions of ion linacs, such as the transformation of radioac- 
tive waste, the production of tritium, and fusion drivers. 
High currents of the order of 100 mA restrict beam losses 
below 1 ppm. Thorough studies are necessary to under- 
stand mechanisms of intense-beam losses, in particular, 
beam instabilities and halo formation. 

Most of the theoretical efforts so far have concentrated 
on the Kapchinsky-Vladimirsky (KV) distribution of par- 
ticles in transverse phase space [|] . The KV beam density 
is uniform so that space-charge forces inside the beam are 
linear. It allows an analytical investigation and results 
are used to predict the behavior of real beams. On the 
other hand, it is recognized that the KV model, in which 
all particles have the same transverse energy, is not a 
realistic beam distribution, e.g. The present paper 
compares the KV beam with other, nonlinear particle- 
density distributions, which can serve as better models 
for real beams. 



II. ANALYTICAL CONSIDERATION 

We study a continuous axisymmetric ion beam in a 
uniform focusing channel, with longitudinal velocity Vz = 
13c. The Hamiltonian of the transverse motion (v^_ <^ Vz) 
is 



H{r, s) = + ky/2 + g$(r)/(m73/32c2) 



(1) 



where m and q are ion mass and charge, fco is the focusing 
strength of the channel, 7 = (1 — /J^)^-*^/^, r = \/ + y'^ 
is the distance from the z-axis in the transverse plane, 
and s = \/ x'"^ + y'^ {x' = x/ [3c, y' = y/(3c) is the dimen- 
sionless transverse velocity. The electric potential $(r) 
must satisfy the Poisson equation 



(2) 



where f{x,y,x',y') = f{r, s) is the distribution function 
in the transverse non-relativistic 4-D phase space. The 
integral on the RHS is theparticle density d{r). 

Since the Hamiltonian (p is an integral of motion, any 
distribution function of the form f{r,s) = f{H{r,s)) is 
a stationary distribution. We consider a specific set of 
stationary distributions for which the beam has a sharp 
edge (for all ions r < Vmax = a), namely. 



Nnn{Ho - i/)"-i for H < Hq 
for H > Ho 



(3) 



The normalization constants Nn are chosen to satisfy 
27r rdrd{r) = /, where d{r) is the particle density, 
and / is the beam current. The set includes the KV 
distribution, /□ = NoS{Hq — H), as a formal limit of 
n — > 0, as well as the waterbag (WB) distribution, 
/i Ni9{Ho — H), where 9{x) is the step-function. For 
a detailed discussion of these two specific examples see 

i- 

We introduce the function G(r) = Hq — H(r, s) + s /2, 



because the density can be expressed from (H) as d(r) = 
27rA^„ G"(r). Physically, this function gives the maximal 
transverse velocity for a given radius, Smax{r) = •\/2G(r), 
and defines the boundary in the phase space (r, s) . It 
allows us to rewrite Eq. (j|) as 



[rG'(r)]'/r- A2G"(r) 



(4) 



with boundary conditions G(a) — 0, and G(0) = Go 
is finite. Here the parameter — K/ [J^ rdrG"^ {r)\ , 
where K = 21/ (IqP^ 7^) is the beam perveance, and 
In = Aireomc^ /q is a constant. Particular solutions to 
d) are easy to find for n = (KV) and n = 1 (WB) [|. 
For n > 2 a numerical solution is required. 

To compare different transverse distributions on a com- 
mon basis, we consider rms-equivalent beams which have 
the same perveance AT, rms radius, and rms emittance £. 
To characterize the space-charge strength, one introduces 
an equivalent (or rms) tune depression 



V = /l^~K/{4l4^) 



(5) 



which reduces to the usual one for the KV beam. For 
numerical simulations we use dimensionless variables: 
z — koz, and x — Xy/ko/£, where £ — a'a. In normalized 
variables the beam matched radius is a = V (Gg /C-2)/r\, 
where Gg = £/£ and G2 = x^^^/c? . For the KV case, 
Gg = G2 = 1/4, so that a = l/^. The "hats" are 
omitted below to simplify notation. 
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III. NUMERICAL SIMULATIONS 

We use particle-in-cell simulations to study beam sta- 
bility and halo formation in the presence of instabilities. 
A leap-frog integration is applied to trace the time evo- 
lution for a given initial phase-space distribution. The 
space-charge radial electric field of an axisymmetric beam 
can be found from Gauss' law by counting the numbers 
of particles in cells of a finite radial grid, which extends 
up to four times the beam matched radius. The initial 
phase-space state is populated randomly but in accor- 
dance with (H) for a chosen n. The matched distributions 
remain stable except for a minor dilution related to nu- 
merical errors. However, even the matched KV beam is 
unstable for rj < 0.4, in agreement with existing theory 
1^ and earher simulations 

The beam breathing oscillations are excited by load- 
ing a mismatched initial distribution = ^fi, r[ — r[/fj-, 
where f^, r[ correspond to the matched one, and the mis- 
match parameter /i < 1. A typical range of the simula- 
tion parameters: time step At ~ T/lOO, where T is the 
period of breathing oscillations, total number of parti- 
cles Npar = 16K to A096K, where K = 1024, and radial 
mesh size Ar = a/128 to a/16. The code performs simu- 
lations of about 100 breathing oscillations per CPU hour 
for Npar = 256K on Sun UltraSparc 1/170. 

The beam behavior is studied as a function of tune 
depression rj and mismatch /z. Due to a discrete filling 
of a mismatched beam distribution in simulations and, 
for n > 1, due to non-linear space-charge forces, higher 
modes are excited in addition to the breathing mode. 
Some of them can be unstable depending on values of rj 
and fj,. A detailed numerical study of stability and halo 
formation for the KV beam and its comparison with the 
theory predictions |^,|l| have already been reported in 
10,^. Here we compare results for different transverse 
distributions. In Figs. 1-2 the maximal radius of the 
whole ensemble of particles is plotted versus the number 
of breathing oscillations for the KV and WB beams, for 
the particular case oirj = 0.7 and /i = 0.8 {Npar ~ 256K, 
Ar = a/64:). 
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FIG. 1. KV beam radius versus the number of breathing 
periods for rj — 0.7 and /j, — 0.8. Stars are for period averages, 
dots show minimum and maximum during a period. 
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FIG. 2. The same as Fig. 1, but for the WB beam. 

Comparison of Figs. 1 and 2 shows that for these pa- 
rameters the WB beam remains stable much longer than 
the KV one, but eventually it also blows up and some 
particles form a halo far from the beam core. Results 
for the n = 2 distribution are similar to those for WB. 
Some results depend on simulation parameters; e.g., it 
takes a smaller number of the breathing periods for a 
beam to blow up if Npar is smaller (i.e., higher noise). 
However, the maximum radius, as well as the fraction of 
particles outside the core, are practically independent of 
Npar- The number of particles which go into the halo 
and produce jumps of Tmax seen in Figs. 1-2, might be 
rather small. We define the halo intensity h as the num- 
ber of particles outside the boundary — 1.75a divided 
by Npar- Such a definition is arbitrary, but convenient 
to compare beam halos over a wide range of tune depres- 
sions. While the beam behavior in Figs. 1 and 2 seems 
qualitatively similar, the halos for these two cases are 
very different: h ~ 3.5 • 10"'^ for KV, and about 100 
times less for the WB, with only a few particles in the 
halo (less than 10 of 256K). That is the reason for os- 
cillations of Tmax in Fig'- 2: these few halo particles can 
initially all come back to the core simultaneously. 
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FIG. 3. Beam behavior (qualitatively) versus tune depres- 
sion rj and mismatch fi. 

A qualitative picture of the beam behavior for various 
values of the tune depression and mismatch is shown in 
Fig. 3, and is practically the same for all distributions 
studied. 'H' corresponds to beam instability with halo 
formation, usually with a noticeable emittance growth, 
'U' means that the beam is unstable but a halo is not 
observed in our simulations, and 'S' indicates beam sta- 
bility. The most surprising feature of the diagram is the 
lack of any significant dependence on rj for mismatched 
beams; on the contrary, the qualitative changes depend 
primarily on /i. When /i changes from 0.6 to 0.8, the 
ratio £fin/£ini decreases from 1.7-2 to 1.03-1.07 for the 
KV beam, and from 1.4-1.5 to 1.00-1.01 for the WB and 
n = 2. The number of breathing periods after which 
the beam radius starts to grow noticeably and the halo 
forms, has some dependence on 77; it is smaller for small 
V- 

We performed a systematic study of the KV, WB, and 
n — 2 distributions for tunes rj from 0.1 to 0.9 and mis- 
matches fi from 0.6 to 1.0. Figure 4 shows the ratios of 
the halo radius to that of the matched beam for the KV 
and WB beams with three different mismatches, fi = 0.6, 
0.7, and 0.8. Results for n = 2 beam are not shown; they 
are slightly lower than those for the WB beam. The 
KV halo has a larger radius, especially with small space 
charge (large ry), but for space-charge dominated beams, 
at very small rj, the ratios converge for all distributions. 
The analytical model for the KV halo formation ||l| pre- 
dicts finite values of r„iax/a- between 2 and 2.5 depending 
on rj and ji. One can see from simulations that it works 
well also for WB and n = 2 beams. 
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FIG. 4. Ratio of halo radius to that of the matched beam 
for KV (top curves) and WB (bottom curves) beams ver- 
sus tune depression r) for different mismatches: = 0.6 
long-dashed, /i = 0.7 short-dashed, /i = 0.8 dotted. 

Simulation results for halo intensity h are shown in 
Fig. 5 for KV and WB distributions. Again, results for 
n = 2 are just slightly lower than for the WB beam, and 
not shown. The intensity depends essentially on the mis- 
match, and decreases quickly as the mismatch decreases. 
The WB halo is about 2-3 times less intense than the 
KV halo for small space charge and large mismatch (0.6 
and 0.7) but, for space-charge dominated beams, the in- 
tensities are about the same. For /i = 0.8, however, the 
WB halo is at least an order of magnitude less intense 
than the KV one; it is not even included in Fig. 5. An 
apparent decrease of ft. as 77 decreases is due to the defini- 
tion used: the halo boundary radius = 1.75a increases 
as 1/y^. If a fixed boundary is used instead, the same 
for all tunes, the halo intensity would be larger for larger 
space charge. 




FIG. 5. Halo intensity for KV (solid) and WB beams 
(dashed) vs tune depression 77 for mismatches /i = 0.6 (top 
pair), /i = 0.7 (middle pair), /i = 0.8 (bottom, KV only). 



One more interesting feature is how fast the halo de- 
velops. For the KV beam, the process is usually rather 
fast, and the halo saturates after a few hundred breath- 
ing periods. For the WB and n — 2 distributions, it con- 
tinues to grow rather slowly, and asymptotic values are 
usually reached after a few thousand breathing oscilla- 
tions; it takes especially long for rj < 0.3. Data plotted in 
Fig. 5 correspond to the asymptotic values, after N=5000 
breathing periods for WB and after N=:600 for KV (ex- 
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cept 77 < 0.2, where KV results are also for N=5000). 
These 5000 breathing oscillations correspond to 5-10 km 
of the length for a typical machine, much longer than 
existing proton linacs. 



IV. CONCLUSIONS 

Our simulations show the qualitative similarity of the 
beam behavior for all transverse distributions studied. 
The KV beam can be considered as an extreme case com- 
pared to the WB and n — 2 distributions which are closer 
to real beams. The halo intensity is a few times higher 
and saturates faster for the KV distribution than for the 
other two. 

An interesting new observation is that for axisymmet- 
ric beams under consideration the beam stability and 
halo formation depend primarily on the mismatch, not 
on the tune shift. The halo was clearly observed only for 
large mismatches, at least 20%, and its radius is in agree- 
ment with the analytical model ||l| for halo formation. 
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